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ABSTRACT 

We have reconstructed the temperatures and surface densities in the quiescent 
accretion disc in HT Cas by performing a Physical Parameter Eclipse Mapping analysis 
of archival UBVR observations. Using a simple hydrogen slab model and demanding a 
smooth, maximally artefact-free reconstruction, we derive a formal distance to HT Cas 
of 207 ± 10 pc, significantly larger than the 133 ± 14 pc we derive from a re-analysis of 
the data in the literature. 

The accretion disc is small (0.3 — 0.4 Rn) and moderately optically thin but 
becomes nearly optically thick near the white dwarf. The temperatures and surface 
densities in the disc range from 9 500 K and 0.013 g cm“^ in the center to about 4 000 K 
and 0.04 gcm“^ at the disc edge. The mass-accretion rate in the disc is roughly con¬ 
stant but - at the derived distance - uncomfortably close to those which would prohibit 
the dwarf nova eruptions. 

We argue that the larger derived distance is probably incorrect but is not produced 
by inaccuracies in our spectral model or optimization method. The discrepancy can 
be resolved if the emission regions on the disc are patchy with a filling factor of about 
40% of the disc’s surface. This solves the problem with the high effective temperatures 
in the disc - reducing them to around 6 500 K within a radius of 0.2 - and reduces 

the derived temperature of the white dwarf and/or boundary layer from 22 600 to 
15500K. 

The viscosity parameters a derived from all reconstructed temperatures and sur¬ 
face densities are of order 10-100 and cannot be lowered significantly by invoking a 
lower distance or the filling factor. This situation is easily explained using the same 
patchy nature of the emitting material, since the quiescent disc cannot consist of opti¬ 
cally thin regions alone, but also of a dark and hence cold and dense disc which could 
easily contain most of the matter. If we require global values of a of order 0.1, the 
implied total surface densities are 1-100 gcm“^ - just like those expected for quiescent 
discs awaiting the next eruption. 

We discuss several possible sources of the chromospheric emission and its patch¬ 
iness, including irradiation of the disc, thermal instabilities, spiral-wave-like global 
structures, and magnetically active regions associated with dynamo action and/or 
Balbus-Hawley instabilities. 

Key words: binaries: eclipsing - novae, cataclysmic variables - accretion, accretion 
discs - stars: HT Cas 


1 INTRODUCTION 

Cataclysmic variables (CVs), are semi-detached binary stars 
consisting of a Roche lobe-filling late-type dwarf (the sec¬ 
ondary star) which loses matter onto a white dwarf (the pri- 
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mary). Unless the primary has a considerable magnetic field, 
the transfered matter has enough specific angular momen¬ 
tum to create a gaseous accretion disc around this accreting 
star. Some CVs - the dwarf novae - have discs which oc¬ 
casionally have luminous phases of high mass-accretion but 
spend most of their time in a faint quiescent state. In qui¬ 
escence, the optical light is usually dominated by optically 
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thin line emission from the disc and continuum emission 
from the white dwarf, the disc and the bright spot caused 
by the impact of the transferred material onto the disc. 

While spectral and eclipse observations of erupting 
dwarf novae show that the disc is optically thick during 
the outburst state, there is considerable uncertainty about 
the state of the quiescent discs. Theoretical models for the 
outbursts invoking a thermal instability as the cause of the 
eruption require that the quiescent disc store up material for 
the next eruption in a less viscous state in which the disc 
should be cool and optically thick (e.g. Ludwig et al. 1994). 
The inner disc may be emptied by irradiation from the hot 
white dwarf (Leach et al. 1999) or by a siphon flow fueled by 
a hot corona in the inner disc (Meyer & Meyer-Hofmeister 
1994). Thus, the expectation is that the outer disc is cool 
and optically thick and the inner disc is warm/hot and op¬ 
tically thin. 

HT Cas was once called “the Rosetta Stone of dwarf 
novae”, because - being one of the relatively few eclipsing 
dwarf novae - it more readily revealed its secrets (Patter¬ 
son 1981). It has seldom been caught in outburst because of 
the unusually long outburst period of about 400 ± 50 days 
(Wenzel 1987). The quiescent periods may even last up to 
almost 9 yearsQ indicating an extremely low mass accretion 
rate in quiescence. The accretion disc is hardly visible dur¬ 
ing quiescence and the bright spot, the impact region of 
the accretion stream from the secondary, is nearly absent 
(Wood, Horne & Vennes et al. 1992, hereafter WHV92). 
The light curves predominantly show the eclipse of the white 
dwarf. The optical spectrum of HT Cas (Young, Schneider & 
Shectman 1981) shows double peaked emission lines of H I, 
He I, He H and Ca H, typical for a high inclination system 
(i = 81°: Horne, Wood & Stiening 1991, hereafter HWS91). 
Consistent with the broad-band photometry, the averaged 
profiles do not show any significant blue/red asymmetry, 
which would normally be attributed to the bright spot. 

WHV92 investigated the accretion disc using multi¬ 
colour eclipse photometry from which the contribution of 
the white dwarf had been subtracted. They produced in¬ 
dependent classical eclipse maps for each colour and then 
modeled them with a two-parameter model spectrum yield¬ 
ing the temperature and surface density distribution in the 
disc. While this analysis suggested that the disc is optically 
thin - as expected - it suffers in detail due to the indi¬ 
vidual smearing of each eclipse map (see Vrielmann, Horne 
& Hessman 1999; hereafter VHH99). Furthermore, the pre¬ 
subtraction of the white dwarf contribution and the use of 
simple Cartesian reconstruction maps can make it difficult 
to produce a consistent physical picture of the various com¬ 
ponents. 

Our present investigation is aimed at improving our 
understanding of the physical state of the quiescent ac¬ 
cretion disc in HT Cas. We apply the Physical Parameter 
Eclipse Mapping method (PPEM) - a distinctly different to¬ 
mographic ansatz (VHH99) which also uses the Maximum- 
Entropy Method (MEM; Skilling & Bryan 1984) - to the 
same UBVR photometry and derive the white dwarf tem- 


t Since it is a circumpolar object for many northern observato¬ 
ries, the likelihood of detecting outbursts is relatively high. 


UBVR Light curves of HT Cas, = 1.75 



Figure 1. The averaged quiescent light curves in UBVR together 
with the fits made simultaneously in all four filters. The contact 
phases of the white dwarf as given by Wood, Horne & Vennes 
(1992) are drawn as vertical dashed lines. The light curves are 
offset by 1.5 mjy, the horizontal dotted lines give the zero level 
for each light curve. In relation to this, the uneclipsed components 
and the residuals are plotted (solid lines). 

perature and the disc structure by simultaneously fitting all 
available eclipse light curves. 


2 THE DATA AND SYSTEM PARAMETERS 

Since the dataset has been previously published and anal¬ 
ysed, details about the acquisition can be found in HWS91. 
The photometric data consist of 16 eclipses taken in Septem¬ 
ber 1982 and November/December 1983 at the Palomar 60- 
inch telescope using the Stiening photometer which gathers 
data in four filters simultaneously. The previous outburst 
occured in March 1980 and the quiescence lasted until the 
super outburst in January 1985. HT Cas was therefore in 
the middle of a quiescence period. 

The original data have a time resolution of 1 s. As can 
be seen in Fig. 1 in HWS91 the basic eclipse shape appears 
to be constant, we do not even see a distiguished difference 
between the 1982 and 1983 data. But the individual eclipses 
suffer heavily from flickering. We averaged all light curves 
in order to enhance any stable patterns and minimize all 
random features, like flickering. The averaged light curves 
were then phase-binned to Aip = 0.0015 near the eclipse 
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(within phases —0.07 and 0.07) and Aip = 0.003 outside 
this range (Fig. ^). 

We used the photometric system parameters as given 
by HWS91, namely a mass ratio q = Mrdj■M.-wd = 0.15, 
an inclination angle i = 81°, as well as a white dwarf mass 
and radius of Afu,d = O.BlAf© and R^d = O.OllST?.© = 
0.0263i?i,i. With an orbital period of 106 min, the distance 
i?Li from the primary to the Li-point is 3.13 x 10^° cm. The 
mass and radius of the secondary star derived by HWS91 
are 0.09 TW© and 0.154 77©. 

WHV92 determined the E{B — V) as dehnately smaller 
than 0.2 and more likely between 0 and 0.1. Therefore, we 
neglect the interstellar reddening for our calculations. 


3 DISTANCE ESTIMATES 

Wood et al. (1995), using UBVR photometry of the white 
dwarf extracted from the eclipse light curve, derived a value 
of 165±10 pc. Marsh (1990) deduced a distance of 140±14 pc 
using the observed Na I and TiO absorption spectrum of the 
secondary and his own calibration of the TiO AA7165, 7665 
band strengths, a surface brightness calibration for M giants 
(the Barnes-Evans relation; Barnes et al. 1978), and the ra¬ 
dial velocity amplitude A 2 . This derivation makes several 
assumptions about the nature of the secondary star. Marsh 
predicted a mid-eclipse infrared H magnitude of 14.9 ± 0.1 
for the secondary but Berriman et al. (1987) measured 
15.7 ± 0.2, indicating that the system may be farther away 
than 200 pc. Two months after the first infrared measure¬ 
ments, Berriman et al. measured a much brighter H mag¬ 
nitude of 14.5, apparently conhrming Marsh’s distance, but 
this difference could have been due to an increase in the 
constribution from the uneclipsed component (e.g. the outer 
regions of the disc). Although Zhang et al. (1986) derived 
an even larger distance of 320 pc from their analysis of BVR 
light curves, this distance is undoubtably made too large by 
their assumption of pure blackbody emission from the disc 
and the lack of U-band data. 

Marsh’s estimate depends on his identification of the 
M5.4 ± 0.03 spectrum (derived from his own calibration 
of the TiO band strengths), the secondary’s flux at mid¬ 
eclipse (0.224 ± 0.01 mjy at 7 500 A), and the Barnes-Evans 
relation (the surface brightness of giant stars versus V-R 
colour); Barnes et al. (1978). Beuermann et al. (1999) have 
recalibrated the Barnes-Evans relations (surface brightness 
versus colour) for both giants and dwarfs and Hnd that the 
relation for M5 dwarfs is about 0.12 mag fainter than that for 
giants. This effect implies a 5% decrease in Marsh’s distance 
estimate (132 ± 13 instead of 140 ± 14 pc). Marsh’s spectral 
type depends upon the relative strengths of the TiO features 
at AA7165, 7665: the latter band is contaminated with weak 
01 emission, so that the A7165/A7665 ratio may be over¬ 
estimated, and the spectral type later than M5.4. Though, 
the type is probably not any later than M6, since the region 
around TiO A7165 does not show the multiple absorbtion 
features prominent for later types. 

Weichhold & Beuermann (2000) and Beuermann (2000) 
have also constructed a TiO surface brightness relation 
which can be used for systems like HT Cas with severe veil¬ 
ing. Using Marsh’s M5.4 db 0.3 spectral type and assuming 
Mrd = 0.09 Mq, the observed TiO flux at mid-eclipse and a 


94% correction for the Roche flattening at phase 0.0 implies 
a distance of 135 ± 14 pc, in agreement with our revision of 
Marsh’s estimate. If the spectral type is as late as M6, this 
number drops to 93 pc, since the flux of the TiO bands drops 
rapidly after M5.5. 

One can reverse the calculation and compute the ex¬ 
pected spectral type given the observed TiO flux and as¬ 
sumed distance of 165 pc (Wood et al. 1995): this results 
in a spectral type of M4.5. This would imply an increase in 
the A7165/A7665 ratio of 43% over that measured by Marsh 
(1990), a ratio which is already too large due to 01 contam¬ 
ination. A spectral type as early as M2-M3 (corresponding 
to a distance of ~200 pc) can be simply ruled out because it 
would require twice as much light at mid-eclipse as is avail¬ 
able for the secondary and the residual disc. The only way 
to change these results - given that the constraints come 
from the unilluminated backside of the secondary star - is 
to decrease the strength of the TiO band by lowering the 
metalicity to at least \M/H] « —1.0 (to get a reduction of 
the band depth by a factor (133/200)^ = 44%). However, 
there is no evidence to date of any low-metalicity CV’s (see 
Beuermann et al. 1998). An earlier spectral type is also im¬ 
plausible due to the fact that there is weak evidence that 
the spectral types of secondary stars in CV’s have later - 
not earlier - spectral types than normal main-sequence stars 
(Beuermann et al. 1998). 

We must conclude that Marsh’s observations clearly 
suggest a distance of 133±14 pc (weighted average of 132±13 
and 135 ± 14pc). This result remains in conflict with the 
eclipse-mapping results of Wood et al. (1995). Given that 
the latter depends upon several assumptions about the vis¬ 
ibility of the white dwarf and its temperature, and that the 
plausible corrections to our re-analysis of Marsh’s result tend 
to make the distance smaller, not larger, we feel that our re¬ 
vision of Marsh’s estimate is still more reliable. Still, given 
the uncertainties, it is possible that the distance is somewhat 
larger. 

In principle, the PPEM method can provide us with an 
independent estimate of the distance, as previously shown 
by Vrielmann, Stiening & Offutt (2002) if a proper spectral 
model is used. The use of an improper distance produces 
characteristic features in the maps of the fitted quantities. 
In Appendix ^ we demonstrate that reliable distances can 
be found by maximizing the entropy of the fit (making the 
disk maps as simple as possible) at fixed (requiring a 
good fit to the lightcurves). When the wrong distance is 
used, the disk maps become more complicated, and the en¬ 
tropy accordingly becomes increasingly negative. The best 
distance is thus identified as that which maximizes the en¬ 
tropy (simplicity) of the disk maps. The determination of 
a self-consistent distance to HT Cas was hence one of the 
objectives of this study. 


4 THE PPEM ANALYSIS 

We used the two-parameter spectral model described in 
VHH99 to map the temperature T and surface density E 
of the disc material. The assumed spectrum is that of a sim¬ 
ple uniform slab, 

7, = H.(r)(l-e“""/°°"') (1) 
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where Tv = pK^ 2H is the optical depth, p = T, j 2H is the 
mass-density, H « Csj^Kepier is the vertical scale height of 
the slab, here assumed to be the disc’s semi-thickness, k is 
the mass absorption coefficient, Cs is the sound speed and 
^Kepler the Keplerian angular velocity of the disc material. 

Although we assume a very simple pure hydrogen slab 
in local thermodynamic equilibrium (LTE) with no line 
emission, the difficulty in computing accretion disc spectra 
and our ignorance of the true vertical structure of optically 
thick or thin discs makes this model still very useful as a 
means of distinguishing between optically thick and thin 
parts of the disc. 

At the densities and temperatures considered here, LTE 
is a very good approximation. The assumption of pure hy¬ 
drogen is quite adequate for hot or warm regions of the 
disc. For realistic slabs (i.e. including metals) at temper¬ 
atures lower than about 6 300K, the opacity is still mostly 
provided by H“ and hydrogen bound-free emission but the 
total opacity is much higher due to the increase in the num¬ 
ber of free electrons relative to a pure-hydrogen gas. Thus, 
our method severely underestimates the intensity emerging 
from cooler parts of the disc. If a certain flux is required, the 
algorithm will be forced to increase the fitted temperature 
(which increases the opacity in this range) and/or increase 
the surface density. Regions of the reconstructed disc with 
kinetic temperatures less than 6 300K will thus - in fact - 
have true temperatures and/or surface densities which are 
significantly lower than those reconstructed using our simple 
model. Fortunately, this effect will only be important in the 
least-constrained outer disc where our ability to reconstruct 
physical parameters is the worst anyway (due to the dis¬ 
appearance of the Balmer Jump and low intensities in the 
observed filters). We will see that our upper limits on the 
parameters in some regions are still physically useful. 

For theoretical spectra with no line emission, the physi¬ 
cal parameter maps derived from the UBVR magnitudes are 
similar no matter whether one convolves the spectra with the 
filter band-passes or uses the monochromatic fluxes at the 
mean wavelengths of the pass-bands (indeed in the latter 
case our numerical results are slightly better). This is due 
to the fact that the continuum is very flat (similar fluxes in 
all four passbands). The inclusion of emission lines would 
improve the sensitivity of the maps to the U-band fluxes 
and regions with substantial Ha, and hence to the optically 
thin parts of the disc. However, this effect goes in the same 
direction as that mentioned before for the continuum: these 
less realistic models have less flux in the filter pass-band and 
so require higher temperatures and surface densities. 

Since HT Cas was in quiescence at the time of observa¬ 
tions, we assumed a geometrically thin disc for our PPEM 
study - an assumption which we will later verify. 

We model the white dwarf emission using theoretical 
white dwarf spectra and treat the star as a true spherical 
star which is occulted by the secondary, occults the accretion 
disc behind itself, and can be preset to be partially occulted 
by the accretion disc. Wood & Horne (1990) find that the 
white dwarf is not likely to be partially occulted, although 
they could not completely rule out a partial occultation of 
the white dwarf. Our reconstructions (Section ^) imply 
that the inner disc is optically thick and therefore would 
lead to an occulted “lower” hemisphere of the white dwarf. 
On the other hand, evaporation of the inner disc (Meyer & 


Meyer-Hofmeister 1994) could lead to a small hole in the ac¬ 
cretion disc which, with our limited spatial resolution, would 
be difficult to detect even with PPEM. However, for the 
white dwarf to be completely visible, the hole must be larger 
than rhoie = i?mdtani = G.SR^d ~ 0.17 Rli, which should 
be easily detectable in the reconstruction. Thus, following 
Wood & Horne (1990), we first chose a fully visible white 
dwarf but then also reconstructed the accretion disc using a 
partially white dwarf for comparison (see Appendix . 

We assume the white dwarf has a uniform tempera¬ 
ture T.ujd, with no limb darkening (u = 0, Wood & Horne). 
Therefore, we use HWS9rs white dwarf radius of 0.0263J?);,!, 
which is derived for a vanishing limb-darkening coefficient. 
Since the mean wavelengths of the filter pass-band partly 
coincide with the absorption lines in the theoretical white 
dwarf spectra, we integrated these spectra over the photo¬ 
metric filter pass-band functions to derive the theoretical 
white dwarf flux. These were then used within PPEM to 
determine the white dwarf temperature. 

The slow ellipsoidal variations of the secondary star dur¬ 
ing quiescence and over the small range of orbital phases 
considered is roughly ~13% of the secondary’s flux (Marsh 
1990): given that the secondary star contributes certainly 
no more than 10% of the light in the R-band and much less 
in the others, we neglect the variation in the contribution 
of the secondary over the limited phase range between -0.1 
and 0.1 and treat the secondary’s light as simply part of the 
fitted uneclipsed component. 

For our reconstructions we used a default map with 
medium blurring, as described in Appendix ^ 


5 RECONSTRUCTIONS USING THE 

D ~ 150 PC ESTIMATES 

We first attempted to obtain PPEM paramater maps of the 
disc in HT Cas using the short distances of d = 140 — 165 pc 
estimated by Wood et al. (1995) and Marsh (1990). 

Unlike the results of Wood et ah, our reconstructions 
must fit all of the data simultaneously and so suffer less from 
additional arteficial spatial smearing. Using d = 165 pc we 
see typical features in the maps arising when too small a 
distance was chosen (Appendix ^ : they show emission in 
two bright regions along an axis perpendicular to the binary 
axis with one of them showing a curiously hot pixel. This is 
the algorithm’s desperate attempt to preserve the width of 
the eclipse profiles while minimizing the emitting area. The 
reconstructed disc usiug d ~ 150 pc would simply produce 
too much emissiou if it were axially symmetric. 

We could fit the data with a reduced /N of 2.5 using 
a distance of 165 pc, but the entropy of the solution is very 
low compared to the entropies of the maps corresponding 
to larger trial distances, i.e. there are more arteficial struc¬ 
tures and poorer fits for d ~ 150 pc. Reconstructions using 
Marsh’s smaller distance are even worse. 

Otherwise, the reconstructions show similar radial tem¬ 
perature, surface density and derived parameter profiles as 
Wood et al. (e.g. a flat 8 000 K temperature profile within 
Q.IIRli except for a nearly optically thick arteficial spot 
with a temperature of 10 000 K) and a disc radius of 0.35i?Li. 
The white dwarf is reconstructed with a temperature of 
about 19 300 K, close to Wood et al.’s 18 700T1 800 K. Note 
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Figure 2. The entropy of the reconstructions (a measure of 
smoothness) as a function of the assumed distance. For each trial 
distance the data were fitted with 'x^/N = 2.5 (triangles), 2.0 
(squares) and 1.75 (circles). The dashed lines are parabolic fits to 
the data, they peak at 214 pc, 210 pc and 207 pc, respectively. 
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that this reconstruction contains disc temperatures which 
should be reliable (i.e. are larger than 6 300K). The effec¬ 
tive temperature is also close to expected values of about 
7 000 K at small radii. The optical depths range between 0.1 
to 0.2 except near the white dwarf, where they reache values 
of ~0.7. 


6 RECONSTRUCTIONS INCLUDING A 

FITTED DISTANCE 

Given the problems in the reconstructions for d ~ 150 pc 
estimates, we next attempted to find the distance to HT Cas 
by using the entropy information in the maps (Vrielmann et 
al. 2002). 

We reconstructed the accretion disc maps for a set of 
trial distances and determined the entropy for a fixed /N. 
The result is presented in Fig.^. A parabolic fit to the func¬ 
tion for yA ~ 1-75 peaks at 207 pc, a distance which is sig¬ 
nificantly larger than the previous estimates using the white 
dwarf or the secondary flux. As a confirmation for this value, 
in Appendix ^ we show the maps using a larger distance of 


240 pc showing typical features of too large a trial distance 
used. 

We estimated the error in the PPEM distance estimate 
by performing the test described in Appendix ^ using an 
arteficial disc with a dominant bright spot. The resulting 
distance is 2.5% larger than the original one. This translates 
into an error of at least ±5 pc for our distance of 207 pc. 

In the following Sections we present the maps and de¬ 
rived quantities using the distance of 205 pc which should 
effectively give the same results as 207 pc. We use this dis¬ 
tance although it does not agree with previous estimates 
and although we do accept the smaller dist ance of 133 pc 
because - as we will explain in Section 7.2 - the maps do 
represent part of the real disc. 


6.1 The light curves 

Despite the simplicity of the hydrogen model, we find a very 
good fit to the eclipse light curve, as shown in Fig. with 
a /N of 1.75. The residuals (at the bottom of the plot) 
show little correlation, apart from the flickering between 
phases 0.055 and 0.075. In the U filter, this looks like the 
re-appearance of a bright spot, but the light curves in the 
filters V and R rather suggest two flares at phases 0.055 and 
0.07. 

The residuals between the observed and fitted light 
curves are at most 5.6% throughout the eclipse (between 
phases -0.05 and 0.05), in general largest for the V light 
curve and smallest for U and B, i.e. the Balmer Jump is 
very well reconstructed. The larger residuals of the V and R 
light curves may be due to the lack of line emission in our 
model, particularly Ha. This is also compatible with the R 
light curve having slightly negative residuals in mid-eclipse. 


6.2 The reconstructed accretion disc 


The reconstructed temperature T and surface density E dis¬ 
tributions (maps) are displayed in Fig. ^ and Fig. |. The 
most prominent features are the radially decreasing tem¬ 
perature and increasing surface density. However, before re¬ 
liable quantitative results can be drawn from these maps, 
the (T, E) parameter space has to be analysed in a way de¬ 
scribed in detail in VHH99. 

This analysis shows that both the temperature and sur¬ 
face density values are very well constrained up to a radius 
of 0.25i?i,i. For radii larger than 0.25i?ii, the fitted tem¬ 
peratures and surface densities become significantly anti¬ 
correlated and hence ambiguous. Since we reconstruct both 
parameters simultaneously, we cannot rely on the exact val¬ 
ues of these two parameters for radii larger than 0.25i?Li 
and any derived parameters (with the exception of the ef¬ 
fective temperature, see Section 3.5). To illustrate the reli¬ 
ability of the reconstructed values of T and E, Fig.^ shows 
azimuthally averaged parameter distributions with approx¬ 
imate error bars calculated using the assumption that the 
spectrum at each particular point may vary by up to a /A 
of 3. This yAIN is larger than that of the fit to the data and 
therefore allows some deviation of the parameters due to 
MEM smearing. 

Given local values of temperature and surface density, 
we calculate the specific intensity distributions and find that 
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Figure 5. The reconstructed temperature (left) and surface density (right) map. The Roche-lobe of the primary component lies just 
about outside the plotted region. The dashed line at 0.25iJr,i indicates that for larger radii the reconstructed parameter values become 
ambiguous, the solid line marks the disc edge. The curved, dotted lines in the grey-scale plot are theoretical accretion stream paths for 
mass ratios 0.15 it 50%. The secondary is at the bottom. 


the inner disc’s R-band intensity is 30% of that of the white 
dwarf. This drops with radius to 7% at 0.25 Rli, and to 
less then 1% at 0.4 Rli which we call the disc edge. The 
region outside 0.25 Rli contributes so little flux that the 
data provide only upper limits: the T and E values found in 
this range are controlled by the entropy extrapolating the 
radial T and E gradients defined by the data for the region 
inside 0.25 Rli- 

The temperature in the inner part of the disc lies at 
about 9 500 K (with a range of 8 500 K to nearly 11 000 K) 
and drops to about 7 000 K at 0.25i?r,i and to 4000 K at the 
edge of the disc. The surface density increases in the same ra¬ 
dial range from about 0.013 gcm“^ to about 0.026 gcm“^ at 
a radius of 0.25i?ii and 0.038 gcm~^ at a radius of 0.4i?i,i. 
Given that the temperatures for R > 0.25i?ii are lower than 
6 300 K, - corresponding to temperatures where our spectral 
model is bound to underestimate both physical parameters 
- the true temperatures and/or surface densities are likely 
to be lower. Thus, we estimate that the outer disc radius is 
0.3-0.4 Rli. 

To obtain a bolometric measure of the optical depth, 
we compare the spectrally averaged intensity distribution 
/(T, E) to the black body intensity distribution Ibb(T) de¬ 
rived from the reconstructed temperature T (Fig.™. The 
ratios of the intensities at each radius show that all of the 
disc has an optical depth of order unity, but the inner parts 
are more nearly optically thick. A more quantitative picture 
is drawn by the spectra of different disc regions as shown 
in Fig. The parameters are selected according to the az- 
imuthally averaged maps (Fig.^. While the Balmer Jump 
is present in the central parts, it disappears for radii larger 
than 0.25Rii, because of the low intensity at the blue end of 
the spectrum. Note that the disc is optically thin at all radii 
in the B-band, is marginally thin in the V-band, is optically 


thick in the inner disc in the U-band, and nearly optically 
thick in the R-band. 

The leading side of the inner disc has higher optical 
depths and so is closer to a black-body radiator. This coin¬ 
cides with the larger temperature and surface density values 
seen in Fig. ^ at small radii. This asymmetry is reproduced 
in the eclipse profile of the accretion disc. To explain the 
asymmetry in the light curve, we seperate the phases in four 
different ranges A,B,C,D with A < (j).md,i, 4'wd,i < B < 0, 
0 < (7 < (ji^d,^ and (/u,d,e < D where (j).uid,i and denote 

the phases of white dwarf ingress and egress, respectively. 
Before the white dwarf is eclipsed, we see a gradual, linearly 
sloped disc ingress (phase range A). While the central star 
is obscured, the disc ingress profile is short and rounded, 
reaching minimum light before phase 0 (B). The following 
egress (C) is again more gradual and after the reappearance 
of the white dwarf we see the short, steep egress of the re¬ 
maining disc (D). Phases A and C with the more gradual 
profiles correspond to the following lune, phases B and D to 
the leading, brighter lune of the accretion disc. Hence, the 
leading lune has a steeper gradient in the parameter distri¬ 
butions than the following lune. A comparison with a much 
more symmetric reconstruction of the disc corresponding to 
a fit to the data with x^/-^ = 2.5 shows that exactly these 
features in the eclipse profile are responsible for the asym¬ 
metric parameter distributions. At this higher X^/N level 
we see correlated residuals especially inside the white dwarf 
eclipse profile. 

For comparison, we reconstructed the accretion disc us¬ 
ing a partially visible white dwarf. The results are described 
in Appendix ^ We note here only that considering the 
PPEM method, the solution with the partially visible white 
dwarf is inferiour to the above presented solution. 
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Averaged l(T,I) and Igg(T) distributions Blackbody Distribution 



0 0.1 0.2 0.3 0.4 
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Figure 6. Left: The radial distributions of the spectrally and azimuthally averaged intensity distributions of I{T, S) and Ibb{T). Right: 
The distribution of the ratio It,'s:IIbb a® a gray-scale plot (see text for explanation). The Roche-lobe of the primary component lies just 
about outside the plotted region. For both plots: The dashed line at 0.25i?r,i indicates that for larger radii the reconstructed parameter 
values become ambiguous, the solid line marks the disc edge. The curved, dotted lines in the grey-scale plot are theoretical accretion 
stream paths for mass ratios 0.15 it 50%. The secondary is at the bottom. 


6.3 The white dwarf 

Our white dwarf fluxes are 0.55 mjy, 0.51 mjy, 0.43 mjy 
and 0.31 mJy for the Alters UBVR, respectively. Except for 
the U flux, our white dwarf fluxes are smaller than WHS’s 
fluxes (their Table 1). The remaining flux in the BVR filters 
must originate in the inner part of the accretion disc. 

The reconstructed temperature of the white dwarf is 
22 600 K. This value is larger than Wood et al.’s (1995) 
value of 18 700 ± 1 800 K determined from the same data, 
because of the larger distance used. If we use a partially 
obscured white dwarf (Appendix we arrive at a white 
dwarf temperature of T^d, 1/4, = 26 700K, even larger than 
Wood et al.’s (1995) value. 

The fact that PPEM favours a fully visible white dwarf 
is problematic, because this means we should expect a sub¬ 
stantial hole in the inner part of the disc of the radial size 
O-ITRli- Such a hole is not evident in our reconstructions. 
However, Wood & Horne’s (1990) analysis of the white dwarf 
ingress and egress profiles also favours a fully visible white 
dwarf or a white dwarf fully covered by a boundary layer, 
but specifically excludes a white dwarf with m = 0 and an 
occulted lower hemisphere. 

In our study we have not specified a location of bound¬ 
ary layer, but expected it to be reconstructed in the inner 
disc. Since we do not see any indication of a boundary layer 
in the disc, its temperature must be mixed into the white 
dwarf temperature, making it in our model indistinguishable 
from the central object. In this scenario might lie the solu¬ 
tion to the problem, since the boundary layer has a different 
geometrical and radiative structure than a white dwarf. 


The white dwarf or boundary layer appears to have dif¬ 
ferent temperatures at different times, as reported by Wood 
et al. (1995): it varies between 13 000 and 20 000 K, lower 
temperatures associated with low states of the accretion 
disc. The white dwarf and/or boundary layer cools down 
during quiescence when the accretion rate is reduced and is 
heated up again during an outburst. Since the HT Cas data 
were taken 2-3 years after the previous outburst, the tem¬ 
peratures derived using d ~ 200 pc would appear to be too 
high compared with those of white dwarfs in similar dwarf 
novae. However, the current mean mass-accretion rate does 
not necessarily reflect the long-term rate, so that high tem¬ 
peratures could principally be due to a previous phase of 
high mass-accretion rates (e.g. Schreiber, Gansicke & Can- 
nizzo 2000). 


6.4 The uneclipsed component 

The unecllpsed component is the fitted flux in each band¬ 
pass which does not appear to be eclipsed and includes the 
contribution of the secondary star. The uneclipsed compo¬ 
nent comprises no more than 8-10% of the total. It increases 
to the UV, indicating a hot uneclipsed source like a chromo¬ 
sphere above the disc (Fig. reaching so far above the disc 
that it is never eclipsed (i.e. a few white dwarf radii). 

Assuming an M5.4 main sequence star (Marsh 1990) 
and using the Kirkpatrick & McCarthy (1994) flux tables, 
we can estimate the contribution of the secondary star to 
the uneclipsed component (Tab. |^. All the uneclipsed light 
in V and R can be accounted for if we assume the smaller 
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Table 1. The reconstructed fluxes in mjy of the uneclipsed component (.Fun), the total mid-eclipse fluxes ^-^d the the total 

out-of-eclipse fluxes Fout compared with the flux of an M5.4 main sequence star at 133 and 207pc (F(m 5.4)). Columns 5 and 6 give the 
mean and standard deviation of the uneclipsed flux for reconstructions with three different default maps (little, medium and much, see 
Appendix Q to illustrate the uncertaincy of the reconstructed value. 


Filter 

Tout 

Tmid 

Fun 

< Fun > 

a{< Fun >) 

Fm5.4(133) 

Fm5.4(207) 

U 

1.13 

0.127 

0.099 

0.106 

0.002 

- 

- 

B 

0.98 

0.104 

0.078 

0.075 

0.003 

- 

- 

V 

1.00 

0.103 

0.018 

0.010 

0.011 

0.021 

0.009 

R 

1.20 

0.231 

0.124 

0.082 

0.028 

0.095 

0.039 


Spectra for different disc regions 



Wavelength [A] 

Figure 7. Reconstructed intensity spectra for disc radii 0.03 (up¬ 
permost), 0.08, 0.15, 0.22, 0.28, 0.35, 0.42 and 0.48 Rli (lower¬ 
most) with overall intensity decreasing monotonically with in¬ 
creasing radius. The parameters T, E used to calculate the spec¬ 
tra are the averaged values from Fig. ^ Additionally, black body 
spectra are shown for comparison at temperatures 9 900 K (upper 
dashed line) and 3 650 K (lower dashed line), corresponding to the 
inner and outer disc temperatures at 0.03 and 0.48 Rli, respec¬ 
tively; the circles indicate the mean wavelengths of the UBVR 
band-passes. 


distance of 133 pc to be the correct one (see Tab. |^). The 
fact that actually the expected V and R fluxes are some¬ 
what higher than the reconstructed one in this case is due 
to the uncertainty in the determination of the uneclipsed 
flux. However, it could also indicate that the secondary is 
more likely of slightly later type, a suggestion mentioned in 
Section |^. Alternatively, the hot uneclipsed source is bound 
to have some V and R flux, a situation easily accounted for 



Wavelength [/A] 

Figure 8. The uneclipsed component (dashed line) in comparison 
to the mideclipse flux (solid line) and the contribution of a M5.4 
dwarf star at a distance of 133 pc (dash-dotted line) and at 207 pc 
(dotted line). The points with error bars give a mean and standard 
deviation of the uneclipsed flux for reconstructions with three 
different default maps {little, medium and much, see Appendix^ 
to illustrate the uncertaincy of the reconstructed value. 


in the solution with the larger PPEM distance in which the 
secondary contributes less light. 


6.5 Derived parameters 

The reconstructed T-E-maps allow us to derive further pa¬ 
rameters like the effective temperature or mass-accretion 
rate distribution, the apparent viscosity, and the disc height 
and mass. 

The effective temperature Teg 

The effective temperature is calculated as the integral 
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Effective Temperature Distribution 



Radius (R^i) 

Figure 9. The radial effective temperature distribution for a dis¬ 
tance of 205 pc (circles) and 133 pc (dots) assuming a covering 
factor of 41%. The underlying dotted lines are theoretical steady 
state temperature distributions for log A4 = 13 to 18, the one for 
jC\ = lO^^gs”^ is drawn solid for reference. The dashed-dotted 
line indicates the critical mass accretion rates Ma according to 
Ludwig et al. (1994). The vertical dashed line at 0.25 -Rli indi¬ 
cates the radius above which the reconstructed parameter values 
become ambiguous. The vertical solid line marks the disc edge. 


of the intensity of the emission over all directions Q, and 
frequencies u 


/^(r,E) du dn = uTeff 


( 2 ) 


with a the Stefan-Boltzmann constant and 1^, according to 

Eq.| 

Fig.H shows the effective temperature distribution cal¬ 
culated directly from our spectral model and the fitted pa¬ 
rameters. Even though the reconstructed parameters T and 
E are somewhat ambiguous for radii larger than 0.25i?i,i, 
the effective temperature Tes is reasonably reliable, since 
I(r, E) must still be reproduced (the light curves are fit¬ 
ted very well) and the optical flux is a large fraction of the 
bolometric flux. 

Within a radius of 0.2f?i,i the Tefr distribution is rela¬ 
tively fiat, with values around 7 500 K or 10 000 K. The dis¬ 
tribution of Teft is distinctly spatially separated, following 
the pattern of the optical depth (with higher Tea’s corre¬ 


sponding to the optically thick and lower Tea’s to optically 
thin regions). For radii larger than 0.2Rli, the disc appears 
to be compatible with a steady state disc with a mass ac¬ 
cretion rate around 3 x 10^®gs“^ (i.e. 5 x lO^^^TWgyr”^). 

Also shown in Fig. ^ is the critical effective temperatures 
as a function of radius predicted by Ludwig et al. (1994). 
This corresponds to the critial mass-accretion rates above 
which we would expect steady accretion in an optically thick 
disc to occur. Within a radius of 0.2i?i,i our derived Tefr lie 
above this line. Since the data were taken in the middle of a 
quiescent period, the effective temperatures should not have 
reached above the critical values. 

The viscosity 

The most important parameter describing the disc is 
the viscosity u, usually parametrizised as !/ = aCaH (where 
Ca is the local sound speed and H the disc scale height) using 
the Shakura & Sunyaev (1973) a-parameter. In order to be 
able to fit the dwarf nova light curves with current accretion 
disc models, we need values of a around 0.01 in quiescence 
(e.g. Smak 1984). 

Using the standard relation between the viscously dis¬ 
sipated and total radiated flux iT = IvdTl, 

2HaP-nK = / F^du = crTefl (3) 

^ Jo 

where P is the (gas) pressure, we can map the viscosity 
parameter a of the disc material. The derived values lie be¬ 
tween 20 and 230, far larger than those expected in any kind 
of turbulent accretion disc and similar to those derived by 
WHV92. 

The dise height and mass 

The equilibrium disc height ratio, H/R = Ca/0,K {ca is 
the local sound speed and fix the Keplerian angular velocity 
of the disc material) is roughly 0.01, i.e. the disc has a very 
small opening angle of a degree or less. The innermost part 
of the disc and the disc towards the edge are even narrower 
(with an opening angle of approximately 074) giving the 
disc a slightly concave shape in the inner region and convex 
shape in the outer part. In any case, the reconstructed disc is 
extremely thin, justifying (at least internally) our thin disc 
approximation. 

Within a radius of OARli the mass of the disc is 5.4 x 
10^® g (2.7 X 10“^"* Ado). This factor is uncertain by up to 
50% due to the ambiguity of the surface density density 
values in the outer parts. A more robust calculation of the 
disc mass for the inner part of the disc (up to a radius of 
0.25i?ii) leads to 1.3 x 10^® g (6.4 x 1O"^®A40). 


7 THE RESOLUTION OF THE DISCREPANCY 

Our PPEM distance of 207 pc is significantly larger than 
the estimates of 140 ± 14 pc by Marsh (1990), 165 ± 10 pc by 
Wood et al. (4995) and our revised estimate of 133 ± 14 pc 
from Section H. The main problem with the reconstructions 
using d ~ 150 pc is the appearance of characteristic distor¬ 
tions due to too much flux from the nearer disc, a problem 
which can easily be remedied by increasing the distance. On 
the other hand, the PPEM reconstructions at d ~ 200 pc 
also have their problems, notably effective temperatures, 
Tefi, which are higher than those expected for quiescent 
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discs. Both solutions suffer from the fact that the values 
of the viscosity paratmers, a are orders of magnitude higher 
than can be physically explained. 

We therefore need to check quite critically the potential 
pit-falls in our own analysis and - if the method can be 
shown to be reliable enough - to explain the behaviors of the 
different reconstructions using an improved physical model 
of the disc. 


7.1 Reliability of the PPEM method 

If our reconstructions would like to have a distance which 
is too high, we must identify which false assumption has 
produced this result. There are several possible explanations 
for the discrepancy. 

• This effect could be a peculiar side-effect of the PPEM 
approach since it did not show up in classical eclipse map¬ 
ping analyses. However, the latter do not use the constraints 
on the colors and spectral shapes of the mapped regions and 
are thus free to yield even very unphysical solutions. The si¬ 
multaneous use of all the information as well as the physical 
constraints posed by the physical (if simple) spectral model 
should make PPEM more sensitive to such effects. Thus, it 
is not surprising that this problem did not arise in WHV92s 
analysis and there is no a priori reason why PPEM should 
not produce reasonable results. 

• Wood et al.’s distance is based on the white dwarf 
colours alone, which were deduced from the original eclipse 
light curve while our analysis includes the information of 
the accretion disc. Eurthermore, the distance Wood et al. 
fitted is not independent from the temperature fitted to the 
white dwarf fluxes. We believe our method is therefore more 
advantageous. 

• As we have explained, our model for the emission from 
the disc does not necessarily describe the true emissivity of 
the disc, both because we have no a priori knowledge of the 
structure of the disc and because of our use of pure hydrogen 
spectra. However, the reason that our algorithm produces a 
larger distance is that the surface brightness of the disc is 
basically set by the colours and the eclipse width, which, 
for a given orbital separation, corresponds to a maximum 
size of the emitting region. Solutions with smaller distances 
have disc fluxes which are too large for the given spectra 
and areas, i.e. the specific intensity of the disc would have 
to be smaller. Lowering the surface density would decrease 
the intensities, but also changes the colours, since this would 
produce a large change in the optically thinnest part of the 
spectrum. The algorithm therefore reacts by reconstructing 
a disc with much more scatter, i.e. a much smaller entropy 
and a more or less elliptically shaped emission region with 
the long axis perpendicular to the binary axis. This is an 
unphysical response to a physically constrained but primar¬ 
ily numerical problem. Our assumption of pure hydrogen 
spectra already underestimates - not overestimates - the in¬ 
tensities from the cooler outer disc. Thus, we expect that 
the use of improved spectral models would not necessarily 
solve this distance problem. 

• The distance could also be affected by our assumption 
of a pure hydrogen spectral model neglecting any line emis¬ 
sion from hydrogen and other elements. While it is formally 
not difficult to include the lines, the theoretical line fluxes 


are more dependent upon the detailed model of the verti¬ 
cal structure of the disc, are principally more sensitive to 
the diffuse outer layers of the emission region which do not 
contribute much to the total optical depth, and suffer from 
direction-dependent opacity effects due to the shear in the 
disc (Horne & Marsh 1986). Only a test with an apropri- 
atly complicated model would clarify this; in a future paper 
we will use more realistic disc spectra calculated with mod¬ 
els from Hubeny (1991). A preliminary analysis is shown in 
Vrielmann, Still & Horne (2001). We previously noted, that 
when hydrogen dominates the continuum opacity and sup¬ 
plies most of the electrons, then metals have little effect on 
the spectrum. Only in cool disc regions will the metals have 
a noticable effect on the opacity. However, these regions have 
a low intensity. We therefore think that our assumption is 
acceptable and that the lines do not play a major role at the 
high densities implied by the reconstructed maps. 

• Another assumption that might influence the distance 
estimate is our neglect of a dramatic vertical stratification 
of the disc. While the effect on the distance determination is 
also difhcult to estimate and a full consideration lies beyond 
the scope of this paper, we consider a scenario in which the 
disc consists of a cool disc with a hot chromosphere. In this 
case, the main emission source is at the base of the chro¬ 
mosphere. Because chromospheres tend to stabilize around 
10,000 K due to the strong dependence of opacity on tem¬ 
perature when hydrogen is ionized (Williams 1981, Tylenda 
1981), our assumption should not be so bad. 

• And finally, Vrielmann (2001) & Vrielmann et al. (2002) 
were able to determine a PPEM distance for UU Aqr and 
V2051 Oph that agree very well with previous estimates. 
Furthermore, they found no discrepancy to previous distance 
estimates for the systems IP Peg and VZ Scl. This gives 
us confidence in our method, though naturally, it must be 
tested on more systems. 

As the method itself cannot clearly be blamed for the 
discrepancy, we must find another solution to this discrep¬ 
ancy tied to the structure of the disc itself. 

7.2 The solution: a patchy disc 

An attractive way to reduce the flux from the disc without 
substantially changing the colours and total disc area is by 
allowing for a covering factor: the disc need not be uniformly 
covered with emitting material if only a fraction of the lo¬ 
cal disc surface can produce the observed optical flux. Since 
the “chromospheric” emission modelled here is also seen in 
the near infrared in other quiescent dwarf novae (Berriman 
et al. 1985), the remaining area must be filled with “dark” 
material which may only show up at even longer infrared 
wavelengths, i.e. dark patches which are invisible in the op¬ 
tical part of the spectrum. The covering factor we require is 
about C = (133/207)^ = 41%. 

While the presence of such a covering factor C < 100% 
is a plausible solution of the distance discrepancy, it is at 
present not possible to establish if and how the covering 
factor may change with disc radius or azimuth. In future 
analyses where the distance to a CV is well known one may 
consider adding a covering factor C as an additional param¬ 
eter to be mapped. This may be especially viable if we use 
eclipses of the line emission. 
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Fortunately, the physical parameters of the patchy disc 
located at the (true) d ^ 150 pc will not be very differ¬ 
ent from those reconstructed at d ~ 200 pc since the total 
subtended area of the emitting regions (as opposed to the 
total disc area) is the same. This preserves all of the im¬ 
portant spectral properties except the effective temperatures, 
Teff. The latter are derived from the bolometric fluxes and 
the total (rather than the emitting) areas; the Tefi’s need 
to be scaled down to match the lower distance by a factor 
of (0.41)^/^ =80%. The resulting mass-accretion rates in the 
maps are thereby lowered (as can be seen in Fig.|^ to values 
below or at least closer to the critical mass-accretion rate. 

By moving HT Cas closer, the htted temperature of the 
white dwarf is also reduced. In Appendix ^ we show a re¬ 
construction with a fixed covering factor C = 41% for a 
distance of 133pc. The resulting white dwarf temperature is 
15 500 K. For 150 pc and C = 52%, we get 17450 K, close 
to Wood et al.’s (1995) value of 18 700 ± 1 800. 

The Q values will also be reduced for such a patchy disc, 
but only down to values of 10-100. However, the presence of 
additional but optically unseen material solves the problems 
associated with the still unreasonably large a values: if most 
of the viscous transport occurs in the unseen material and 
the small amount of observed “chromospheric” material h- 
nally radiates away this dissipated energy, then reasonable 
values of the “turbulent” viscosity can be invoked. The rel¬ 
ative amount of unseen matter is then at least equal to the 
inverse of the derived a, i.e., of the order of tens or hundreds 
of grams per cm^. As theoretical calculations (e.g. Ludwig et 
al. 1994) show, the hysteresis curve in the effective temper¬ 
ature - surface density diagram used to explain dwarf nova 
outbursts lies just at these surface density values. 

Thus, if at optical wavelengths we only see the chromo¬ 
sphere, we underestimate the surface densities severely and 
therefore overestimate the a’s. A simple estimate shows that 
just such surface densities of 10-100 g cm“^ in an underlying 
cool disc lead to the desired values of a < 1. 

An analysis of single eclipses (see Fig. 1 in HWS91) 
could in principle reveal the pattern, however, the pattern 
would be smeared out due to the differential rotation of the 
disc. Another problem one encounters with mapping single 
eclipses is the flickering, easily hiding just such small steps 
expected in the eclipse profile. The flickering introduces fea¬ 
tures in the map that cannot be distinguished from real 
patterns of a patchy disc. 

Our analysis therefore reveals the disc as an average 
between the bright and dark patches. If the dark patches 
are so cool that they are basically undetectable at the op¬ 
tical wavelength at which the data were taken, the spectra 
(colours) in the disc are determined by the bright patches. 
A time average, however, is lower in intensity than that of 
an active region, due to the averaging over times when the 
patch is present and absent. This leads to the increase in 
distance that we derive using PPEM. 


8 THE SOURCE OF THE PATCHINESS 

The PPEM reconstructions of the intensity distribution of 
the disc in HT Cas imply that the optical emission does not 
cover the entire face of the disc. If proven, this fact would 
place severe constraints on the implied structure of the disc. 


In this model, most matter is carried by an underlying op¬ 
tically thick, cool accretion disc which generates the energy 
finally radiated away in the chromosphere. This underlying 
disc must have surface densities of the order of tens to hun¬ 
dreds of gcm“^ in quiescence in order to explain the dwarf 
nova eruptions (e.g. Ludwig et al. 1994). 

The presence of patchiness would argue against models 
for chromospheric emission which is simply produced either 
by irradiation-induced coronae (Smak 1989) or thermally 
unstable viscous atmospheres (Adam et al. 1984), since these 
processes should produce a uniform chromospheric layer over 
the entire geometrically thin disc. This argument could be 
weakened if the disc were able to maintain a vertically “cor- 
regated” structure like that suggested in systems with spi¬ 
ral waves (Steeghs, Harlaftis & Horne 1997), in which case 
either irradiation or locally viscous and thermal processes 
could produce the required mottled emission. For example, 
the opening angle of spiral waves in discs is supposed to de¬ 
crease with decreasing azimuthal Mach number, so quiescent 
discs may have very tightly wound-up spirals which either 
produce local dissipation or raise material high enough to 
be irradiated by the central disk/star. 

An alternative model is that the patches are similar to 
active regions on the surface of the sun and could be linked 
to the flickering in the disc. These regions could be due to 
magnetic activity in localized regions on the surface or in 
the upper layers of the accretion disc. For example, magnetic 
flux created by dynamo action and/or the Balbous-Hawley 
instabilities driving the viscosity would rise out of the cool 
midplane regions and dissipate most of the viscously gen¬ 
erated energy via magnetic reconnection or similar coronal 
processes. This model would also explain why the emission 
line intensities are proportional to the local orbital frequency 
(Horne & Saar 1991). 

The previous models generally assume that the optical 
emission region is an upper layer above an optically invis¬ 
ible disc containing most of the mass, even locally along a 
vertical cut through the disc. A very different alternative is 
suggested by the observation of “mirror” line eclipses in the 
Paschen lines in the dwarf nova IP Peg by Littlefair et al. 
(2001) at orbital phases around 0.5. Such behavior can only 
be explained by having a disc which is optically thin enough 
to let the light of the secondary star shine partially through 
the disc at phase 0.5. This observation is inconsistent with 
a uniform, optically thick cold disc with surface patches of 
chromospheric emission on top (and bottom) and suggests 
(1) that the patches implied by the PPEM reconstructions 
are exactly those areas of the disc which are vertically truly 
optically thin and (2) that the dense, cold material needed 
to explain both the surface densities necessary to explain 
the dwarf nova eruptions and the anomalous derived vis¬ 
cosity parameters is laterally rather than vertically distinct 
from the chromospheric patches. The fact that the “mir¬ 
ror” eclipses are seen at large Keplerian velocities argues 
against patchiness which is simply due to the “shredding” 
of the outer disc by the accretion stream (weak in HT Cas) 
and the secondary’s tidal forces. Thus, there would have to 
be some thermal or viscous instability which produces the 
mottled structure throughout the disc. 
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9 LIFETIME AND SIZE OF PATCHES 

If the patches are associated with reconnection events 
(disruption of structures) due to MHD turbulence, the 
timescales for the creation of structures will be linked to 
the local Keplerian timescales (tie(r)). 

tK{r) = 2Ti/Q.K{r) 



For HT Cas, tKir) ~ lOsec near the white dwarf and ~ 
17min (< Porb) at the disc edge (0.4f?i,i). 

The timescales for reconnection events can be estimated 
using the timescales for flickering - between a few tens of sec¬ 
onds and several minutes (Bruch 1992), i.e. covering the full 
range of Keplerian timescales in the disc. Flickering could 
either be caused by temperature fluctuations on the Kepler 
timescale (Welsh, Wood, Horne 1996) or caused by a self- 
organized criticality (or avalanche) effects as described by 
Yonehara, Mineshige & Welsh (1997). 

The maximal lifetime of a magnetic structure is deter¬ 
mined by the differential rotation of the disc material. For 
example, whenever a magnetic loop is rooted at different 
radii, r and r -I- Ar, shearing in the disc will tend to destroy 
its coherence on a timescale 

tco « (tK{r)~^ - tKp + Ar)~^) ^ 



tco is of the order of a few minutes in the inner disc and 
several hours (i.e. several hours) at the disc edge. 

The size of the structures is linked to the disc thickness 
at each radius. If the vertical and radial velocities in the disc 
are comparable, then the radial size may be comparable to 
the disk thickness. The azimuth size of the structure will 
be somewhat larger due to differential rotation shearing out 
the structure. 

It will probably be very difficult to map the structures 
with photometric light curves, because single eclipses will 
contain the projection of many structures. These cannot be 
resolved 2-dimensionally with eclipse durations 3> the flick¬ 
ering timescale, because the flickering prevents us from see¬ 
ing the signatures of the structures in the eclipse profile. 
Furthermore, the patches are moving on a Kelper time scale 
generally (much) smaller than the eclipse duration which 
leads to smeared out patches. Averaging over several cycles 
will average over times when the structures are at different 
places (see also Section 7.2). Thus, we can only expect to 
see the structures in fast and high resolution spectroscopy. 
The patches will lead to Hne-scale wrinkles that will vary 
with the local Kepler velocity along the emission line pro¬ 
files (Hoffmann, Hessman & Reinsch 2002). 


10 SUMMARY 

We have applied the Physical Parameter Eclipse Mapping 
(PPEM) method to archival UBVR photometry of the 
eclipses in the quiescent dwarf nova HT Cas in order to 


map the optically thin matter responsible for the optical 
continuum and emission line spectrum. 

Using a spectroscopically determined distance, we were 
unable to obtain reasonable maps of the temperature and 
surface density in the disc: the reconstructions showed char¬ 
acteristic artefacts showing that the surface brightness of the 
disc was too high. By insisting that the disc be maximally 
structureless, we were able to find reasonable solutions at 
a distance of 207pc. The resulting accretion disc has a ra¬ 
dius of 0.3-0.4i?ii = 9 — 13 X 10® cm, kinetic temperatures 
between 4000 K and 10^ K, surface densities between 0.013 
and 0.04 gcm“®, and t < 1 everywhere except near the hot 
central white dwarf (22 600 K). The effective temperature 
is flat in the inner part of the disc and may follow a steady 
state solution towards the disc edge. We derive apparent val¬ 
ues of the viscosity parameter a from the bolometric fluxes, 
temperatures, and surface densities which are orders of mag¬ 
nitude higher than they should be. 

Acknowledging that the distance to HT Cas is unlikely 
to be as high as suggested by our PPEM entropy analysis, 
we explain why the problem is unlikely to be due to a poor 
treatment of the white dwarf or the wrong choice of a spec¬ 
tral model. We are able to rectify both the poor reconstruc¬ 
tions of the models using d ~ 150 pc and the problems with 
those using d ~ 200 pc by invoking emission regions which 
are not distributed uniformly across the face of the disc but 
are patches with a covering factor of only 41%. Since the 
angular extents of the emitting regions in the distant non- 
patchy and close by patchy discs can be identical, the recon¬ 
structed physical parameters of the d ~ 200 pc solutions are 
still valid at d ^ 150 pc with the exception of those quan¬ 
tities which depend upon the bolometric fluxes: the mass- 
accretion rates (or effective temperatures) need to be scaled 
down by 80%. This reduces the local mass-accretion rates to 
just at or below the critical rates for maintaining the dwarf 
nova eruptions for all but the innermost disc. 

The empirical values of a (cf. - ranging from about 10 
to 1000 - for both d ~ 150 pc and d ~ 200 pc are still orders 
of magnitude larger than is physically plausible. Fortunately, 
the patchy chromosphere on the disc surface requires a sub¬ 
stantial mass of very cool material which powers the viscous 
energy finally radiated away in the chromosphere and which 
have surface densities of the order of tens or hundreds of 
gcm“® - just those needed by the disc instability models. 

The patchiness excludes models for disc coronae due to 
local thermal instabilities or to irradiation. If we interpret 
the patchy chromospheric emission regions as magnetically 
active regions, the fact that the dissipation is proportional to 
the angular velocities could be a sign that what we are see¬ 
ing magnetic flux created by magnetohydrodynamical insta¬ 
bilities and/or disc dynamos responsible for the anomalous 
viscosity in discs dissipated in magnetically active regions. 
However, the presence of “mirror eclipses” in another dwarf 
nova implies that the warm chromospheric patches may be 
situated laterally adjacent to rather than simply above the 
cold, more massive regions. 

HT Cas may really deserve its name as “Rosetta Stone” 
of dwarf novae, since it could be the key to understanding the 
anomalous viscosity in accretion discs as due to hydromag- 
netic turbulence. High time and spectral resolution spec¬ 
troscopy using a 10-m class telescope and/or multi-colour 
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observations in a yet to be detected low state could provide 
the final proof for our hypothesis. 

The Physical Parameter Eclipse Mapping proves to be 
a powerful tool to investigate the physics of accretion discs 
- even when it at first appears to produce a wrong result! 
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APPENDIX A: TEST FOR DISTANCE 
ESTIMATE 


In order to test the reliability of the PPEM distance esti¬ 
mate, we performed a test in which we calculated a light 
curve from an arteficial disc and attempted to reconstruct 
it using various different trial distances. 

Our test disc consists of an axisymmetric temperature 
and surface density distribution with a bright, gaussian spot 
at a distance of O.SRli, an azimuth of 20° and a sigma of 
0.02i?i:,i. We added this non-axisymmetric feature in order 
to test if this will influence the distance estimate. 

The artifical light curves were calculated using the same 
phase resolution as for the HT Cas data displayed in Fig.|^ 
i.e. 0.003 outside the eclipse (phases ±0.065) and 0.0015 dur¬ 
ing eclipse yielding 112 data points between phases ±0.1. We 
added noise with a signal-to-noise level of 100. The correct 
distance was assumed to be 205 pc. 

Fig. Al shows the resulting entropy - distance relation. 
A parabolic fit to the data points peaks at 210 pc, i.e. at a 
distance 5 pc larger than the true distance. 

In several applicaton we noticied that when fitting light 
curves, the distance decreases with decreasing (see e.g. 
Fig. ^). This effect is usually only small, about 5 pc for a 
drop in of 1. In applications to real data, one usually en¬ 
counters a lower limit of x^ below which no good fit to the 
data can be found (or the reconstructions start showing arte¬ 
facts). This limit is determined by flickering in the data or 
by the adequacy of the spectral model used (for a discussion 
on the appropriateness of our model see Section resolution). 
We expect therefore, that the true distance will be slightly 
smaller than the derived distance. 


APPENDIX B: RECONSTRUCTION WITH 
SMALLER DISTANCE 

To illustrate the difference between the maps for a trial dis¬ 
tance of 165 pc as determined by Wood et al. (1995) and our 
PPEM distance estimate we plot the ratio of the parameter 
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distance (Rli) distance (R^^) 

Figure A2. The ratio map in temperatures (left) and surface densities (right) for trial distances 165 pc (Wood et al. 1995) and 205 pc 
(PPEM, see Section ^). The Roche-lobe of the primary component lies just about outside the plotted region. The dashed line at 0.25Rli 
and the solid line at OARli are drawn for comparison with Fig.^. The curved, dotted lines in the grey-scale plot are theoretical accretion 
stream paths for mass ratios 0.15 di 50%. The secondary is at the bottom. 



1 

o 


Test 

cose 



X 







- 



- 



- 



- 


o 

- 


o 

- 


X 

- 



- 


LO 

- 



- 



- 



- 




/ 

\ 




- 

O / 

\ 

- 


1 


/ 

^ o 



o 



\ 





o ' 



X 



\ 



CD 

- 

/ 


- 

Q_ 






p 


■ 

/ o 

\ 

■ 



- 

/ 

Q 

- 


o 

- 

/ 

/O 

\ 

\ 

- 


X 


/ 

\ 



a- 

- 


\ 




- 

o/ 

\ 

\ 

- 



- 

/ 


- 


o 

- 



- 


X 

- 


o 

- 


00 

- 



- 



- 



- 


o 

- 



- 


^— 


_^^_ 

_^^_ 



ot180 200 220 240 

Distance [pc] 


Figure Al. The entropy - distance relation for the test case. 
For each trial distance the data were fitted with = 1- The 

dashed line is a parabolic fit to the data, peaking at 210 pc. 


maps in Fig. A2. To allow comparison, we used the recon¬ 
structions for a of 2.5. An asymmetry can clearly be seen. 
Since the map for 205 pc is much smoother than the one for 
165 pc, the asymmetry is mainly due to the latter map. 

The temperature ratio map shows maxima in two areas 
at azimuths ±90°, while the surface density ratio map shows 
a minimum at azimuth 0°. This shows that the algorithm 
tries to maintain the eclipse width and at the same time 
to reduce the emitting area, leading to a shortening of the 
disc along the binary axis and thus to an elongated emission 
structure perpendicular to the binary axis. Such an emission 
structure is usual for a reconstruction for which the trial 
distance is too small. 

Additionally, the reconstruction shows a single hot pixel 
at approximately azimuth —95° and radius 0.15i?i:,i. Such 
hot pixels are usually artefacts, since it is physically unlikely 
to maintain such a sharp gradient in temperature and in such 
a small area. It is especially unlikely since the light curves are 
averaged. Such a hot pixel is usually an additional indication 
for a false trial distance. 


APPENDIX C: RECONSTRUCTION WITH 
LARGER DISTANCE 


For completeness we show the difference between the recon¬ 
structions for trial dista nces 240 and 205pc as a ratio of the 

Again we used a of 2.5 for 


Cl 


parameter maps in Fig. 
comparison. These maps show also an asymmetry, however, 
a distinctly different pattern. Since the map for 205pc is very 
smooth for a x^ of 2.5, the asymmetry is mainly due to the 
map for 240pc. 

The maps show maxima along an axis parallel to the 
binary axis. This can be understood if one considers the 
following: the eclipse width is defined by the eclipse profile 
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distance (Rl-]) distance (R^i) 


Figure Cl. The ratio map in temperatures (left) and surface densities (right) for trial distances 240 pc (Wood et al. 1995) and 205 pc 
(PPEM, see Section ^). The Roche-lobe of the primary component lies just outside the plotted region. The dashed line at 0.25i?r,i and 
the solid line at 0.4Rr,i are drawn for comparison with Fig. ^ The curved, dotted lines in the grey-scale plot are theoretical accretion 
stream paths for mass ratios 0.15 it 50%. The secondary is at the bottom. 


and determines the extension of the accretion disc perpen¬ 
dicular to the binary axis. However, if the distance is too 
large, the disc needs to cover more area than available if the 
disc were axisymmetric. Therefore, the emission region in 
the disc extends along the binary axis without compromis¬ 
ing the eclipse width. This pattern is thus typical for a map 
that was derived with too large a trial distance. 

The reason for the enhanced structure in the disc is the 
same as for the curious hot pixel in the map reconstructed 
with the small distance. 


APPENDIX D: RECONSTRUCTION WITH 
PARTIALLY VISIBLE WHITE DWARF 

Since the inner part of the disc appears to be optically thick, 
the “lower” hemisphere of the white dwarf should be (at 
least partially) occulted. Therefore, we reconstructed the 
disc and white dwarf temperature with presetting the oc- 
cultation of the central object. 

These reconstructions show significantly more scatter 
than the ones shown in Fig.^ for the fully visible white 
dwarf. This scatter is expressed by a much lower entropy 
of these maps compared to the ones using a fully viaible 
white dwarf: the value of the entropy is —1.1 x 10“®, a fac¬ 
tor of 3.5 below the one for the fully visible white dwarf and 
outside the plotted region in Fig.^ Therefore, the PPEM re¬ 
construction favors the fully visible white dwarf, like Horne 
& Wood^990). 

Fig.^ shows the ratio of the two reconstructions. As 
expected, the white dwarf temperature in the case of the 
white dwarf with a smaller visible area is higher {T^d,i /4 = 
26 700K). In addition to showing more scatter, the values 
of the disc temperatures and surface densities are gener¬ 
ally higher. This is due to a combined attempt to reproduce 


sufficient total flux and the correct colours, since a slightly 
hotter white dwarf shows a slightly bluer spectrum. 


APPENDIX E: RECONSTRUCTION WITH 
VARIOUS DEFAULT MAPS 


During the iteration of an eclipse mapping method like 
PPEM, the intermediate map is compared to a default map 
which in turn is adjusted at certain, predefined intervalls. 
The general use of a default map is explained in more detail 
in VHH99. In short: a default map is a smeered out version 
of an intermediate map calculated using certain blurring pa¬ 
rameters in radial and azimuthal direction. 

The default map therefore introduces a certain degree 
of symmetry or rather smoothness. The degree of smooth¬ 
ness depends on the value of the blurring parameter. It is 
important to note that at the usual blurring parameters 
(see below) asymmetries are not completely or even nearly 
completely suppressed. It only leads to a certain degree of 
smoothness in the disc. This is desireable in order to avoid 
artefacts to dominate the map and to choose one out of the 
many possible solutions. 

In order to estimate the influence of the default map 
on the PPEM solution we used a set of default smearing 
parameters as given in Tab. El, Usually, we use a solid arc 
smearing, i.e. we use a constant arc length. This leads to 
enhanced blurring in the central regions. We use this type 
of blurring, because we expect more large scale structures 
in the outer regions (e.g. non-circular disc, bright spot). In 
one case we used solid angle blurring, leading to enhanced 
blurring in the outer regions of the disc. 

The reconstructed temperature (T) and surface density 
(E) maps all show similar absolut values for T and S only 
with a larger scatter in the maps corresponding to the larger 
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Radius (R^^) 


Figure Dl. The ratio of the temperature and surface density reconstructions using a partially visible white dwarf ^wd/ 4 ) a 

fully visible white dwarf (T^ci /25 ^wd/ 2 )- The dashed line at 0.25 -Rli and the solid line at OARli are drawn for comparison with Fig.^. 


Table El. The blurring parameters. All but one default maps are 
calculated with solid arc blurring. For the only exception where 
we applied solid angle blurring, we used the medium blurring 
parameters. 


description 

Angle (°) 

Radius (Rli) 

tiny 

1 

0.04 

very little 

2.5 

0.1 

little 

5 

0.2 

medium 

10 

0.4 

much 

20 

0.8 

huge 

360 

0.4 


blurring parameter values. As an indicator for the random 
scatter, Fig. E shows the entropy for the solutions using the 
blurring parameters as given in Tab. E 

As surprizing as it seems to find a higher entropy for 
smaller blurring parameters, we have a simple explanation 
for this. Although the maps for smaller blurring parameters 
seem to have more scatter (Fig, p^), this is not random. 
The image in the map rather looks more focussed, tending 
to separate the pixels in the map into either hot or cool 
ones. Considering only the cool pixels, the disc is then very 
smooth. This way, artefacts are also enhanced and can easily 
be misinterpreted. We therefore prefer using intermediate 
blurring parameters and have used the ones called medium 
for our reconstructions. 


APPENDIX F: RECONSTRUCTION WITH 
COVER FACTOR C < 1 


For completeness, we reconstructed the accretion disc for a 
distance of 133 pc implementing a covering factor C = 41% 
for the disc emission. The resulting effective temperature 
(Teft) map is shown in Fig. E The Teff’s range mostly be¬ 
low or around the critical limit given by Ludwi^t al. (199 4), 
similar to the estimate for d = 133 pc in Fig. H. Fig. F2 re¬ 
veals the structures in the disc. They are distinctly different 
to those in maps with distances 7 ^ 205pc or constructed 
with small blurring parameters. We do not know whether 
these structures are real or new artefacts introduced into 
the maps. However, it is likely that the covering factor is 
not constant as assumed in this calculation and that a re¬ 
construction with U as a parameter map will result in more 
realistic solutions. 

The white dwarf temperature is reconstructed to T^d = 
ITdSOA, compatible with Wood et al.’s (1995) estimate of 
18 700 ± 1800 K, especially taking into account that they 
used a slightly larger distance. 
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distance (Rl^) distance (R^i) 


Figure E2. The intensity ratio maps for tiny and huge blurring parameters. The Roche-lobe of the primary component lies just about 
outside the plotted region. The dashed line at 0.25Rli and the solid line at 0.4 _R_li are drawn for comparison with Fig.^. The curved, 
dotted lines in the grey-scale plot are theoretical accretion stream paths for mass ratios 0.15 di 50%. The secondary is at the bottom. 
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Figure El. The entropy for reconstructions with various blurring 
para met ers. The plot gives only the degree of blurring in azimuth, 
Tab. ill also gives the amount of blurring in radius. The circles 
connected by the dashed lines are for solid arc smearing, the single 
point corresponds to solid angle smearing. The entropy for huge 
360° blurring is —7.7 x 10“"^ (it is omitted from the plot for 
clarity). 


© 0000 RAS, MNRAS 000, 000-000 















The Patchy Accretion Disc in HT Cassiopeiae 


19 


Effective Temperature Distribution 



Figure FI. The effective temperature reconstructed for a dis¬ 
tance of 133 pc and a covering factor C = 0.41. The underlying 
dotted lines are theoretical steady state temperature distributions 
for logjCi = 13 to 18, the one for A4 = lO^^gs”^ is drawn solid 
for reference. The dashed-dotted line indicates the critical mass 
accretion rates Ma according to Ludwig et al. (1994). The dashed 
line at 0.25Rli and the solid line at OARli are drawn for com¬ 
parison with Fig.^. 


'r-z/'Biackbody Distribution 



Figure F2. The intensity ratio maps for a distance of 133 pc 
and a covering factor C = 0.41. The Roche-lobe of the primary 
component lies just about outside the plotted region. The dashed 
line at 0.25 -Rli and the solid line at OARli are drawn for com¬ 
parison with the plot in Fig.^ (right). The curved, dotted lines 
in the grey-scale plot are theoretical accretion stream paths for 
mass ratios 0.15 di 50%. The secondary is at the bottom. 
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